Single-shot measurement of triplet-singlet relaxation in a Si/SiGe double quantum dot 
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We investigate the lifetime of two-electron spin states in a few-electron Si/SiGe double dot. At the 
transition between the (1,1) and (0,2) charge occupations, Pauli spin blockade provides a readout 
mechanism for the spin state. We use the statistics of repeated single-shot measurements to extract 
the lifetimes of multiple states simultaneously. At zero magnetic field, we find that all three triplet 
states have equal lifetimes, as expected, and this time is ~ 10 ms. At non-zero field, the Tq lifetime 
is unchanged, whereas the T_ lifetime increases monotonically with field, reaching 3 seconds at 1 T. 

PACS numbers: 73.63.Kv, 85.35. Gv, 73.21. La, 73.23.Hk 



The lifetimes of single electron spins in silicon have 
recently been measured to be as long as seconds in Si 
nanodevices, including gated quantum dots and donors 
[IHl], a promising step towards silicon spin qubits. Two- 
electron singlet-triplet states in a double dot can also 
be used as qubits ISHZ], with the advantages that gat- 
ing operations can be fast and that readout depends on 
the singlet-triplet energy splitting, which can be much 
larger than the single spin Zeeman energy at low mag- 
netic fields. The lifetimes of singlet and triplet states 
have been measured in GaAs double dots and were found 
to depend on magnetic field, falling to < 30 /iS at zero 
field [HI HI- In silicon, neither single-shot readout of the 
singlet-triplet qubit states, nor measurement of their life- 
times has been achieved up until now. 

Here we report measurements of the lifetimes of sin- 
glet and triplet states in a Si/SiGe double quantum dot 
at magnetic fields from 1 T to T obtained using single- 
shot read-out. Using pulsed gate voltages, we repeatedly 
alternate the charge detuning so that it first favors the 
(1,1) charge state (one electron in each dot) and then the 
(0,2) charge state (two electrons in one of the dots.) Be- 
cause of Pauli spin blockade, charge transitions to (0,2) 
will only occur when the spin state is a singlet. We per- 
form hundreds of thousands of such cycles and measure 
the presence or absence of charge transitions using real- 
time charge sensing. By analyzing the statistics of such 
data, we characterize multiple relaxation processes simul- 
taneously, in contrast to time-averaged measurements, 
which are only sensitive to the rate-limiting process. At 
zero magnetic field the triplet and singlet state lifetimes 
are between 5 and 25 ms, lifetimes that exceed those 
measured in GaAs by over two orders of magnitude. As 
magnetic field increases, the lifetime of the Tq remains 
essentially constant, whereas the lifetime of the T_ in- 
creases dramatically, reaching 3 seconds at = 1 T. 
These long times are expected because of the small hy- 
perfine coupling and spin-orbit interaction in Si quantum 
dots. 

The device is fabricated on a phosphorus-doped 
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FIG. 1; (a) SEM image of a device identical to the one used. 
Quantum dots are formed at the approximate locations of the 
two circles. Charge sensing is performed by monitoring the 
current Iqpc through a nearby point-contact, (b) Charge sta- 
bility diagram of the double dot showing the detuning voltage 
Vf . (c) Energies of two-electron states as a function of detun- 
ing energy e. T+, Tq and T_ are the (1,1) triplets; the (0,2) 
triplets are higher in energy. The (1,1) and (0,2) singlets Sn 
and So2 are coupled by spin-preserving, inter-dot tunneling. 
A magnetic field separates the triplet energies by Ez — q^bB. 
(A) Time-averaged occupation of the (0,2) charge state P02 at 
i3| I = with 5 kHz square pulses applied along Vj : pulse am- 
plitudes on gates 2 and 4 are related by Vip = — 0.35V2p. The 
pulses drive (l,l)-(0,2) transitions within the dotted triangle. 
The suppression of P02 above the dashed line shows where 
(1,1) to (0,2) tunneling is suppressed by spin blockade. 

Si/Sio.rGco.a heterostructure with a strained Si quantum 
well approximately 75 nm below the surface. Palladium 
surface gates labelled 1-9 in Fig. [T]^a) are used to form 
the double-dot confinement potential (TU]. A thick RF 
antenna (Ti/Au, 5nm/305nm) is also present near the 
dot gates, but is unused in this experiment. All gates are 
connected to room temperature voltage sources via cold 



RC filters, which are at the measurement base temper- 
ature of w 15 mK. Gates 2 and 4 are also AC coupled 
to coaxial lines, allowing them to be pulsed at frequen- 
cies between 100 Hz and 1 GHz. There is an attenua- 
tion of « 50 dB between each gate and the pulse source 
(a Tektronix AFG3252.) Current through the device is 
measured with a room-temperature current preamplifier 
with a bandwidth « 1 kHz. 

Figure [l|b) shows a charge stability diagram in which 
the absolute occupation of the dots was found by emp- 
tying both dots and then counting electrons back in. 
Fig. [T|c) shows the predicted energies of the two-electron 
states near the (l,l)-(0,2) transition as a function of de- 
tuning energy e, where the transition is at e = [IT] . The 
detuning energy is controlled by varying the voltages on 
gates 2 and 4 along , shown in Fig. [l|b) . The inter-dot 
tunnel coupling tc was measured by determining where 
the Sii and T_ states cross at finite B\\. This is shown as 
^mix in Fig.jljc), and depends on both and the curva- 
ture of the avoided singlet crossing. Using this approach 
[B], we find = 2.8±0.3/xeV (677 ± 73 MHz.) 

To measure the spin of a (1,1) state we pulse the sys- 
tem into a spin blockaded configuration P^14j . where 
the ground state of the system is S02 and the (0,2) triplet 
states are higher in energy than all of the (1,1) triplets: 
T_, Tq and T_|_. We characterize the parameters needed 
to reach this configuration by detecting spin blockade 
in the time-averaged measurement shown in Fig. [ijd). 
Square pulses at 5 kHz are applied along V^. The color 
scale in Fig.[T](d) shows the time-averaged probability P02 
of finding the system in (0,2) as a function of pulse am- 
plitude and offset along V^. When the pulse crosses the 
(l,l)-(0,2) transition, tunneling between charge states re- 
sults in < P02 < 1. The region where this occurs 
is bounded by the dotted triangle in Fig. [l](d). Spin 
blockade occurs in the part of the pulse triangle that 
is above the dashed white line in Fig. [l](d) . Here we see 
< P02 < 0.5, because the system is residing in (1,1) the 
majority of the time. 

Spin blockade does not occur below the white dashed 
line in Fig.jljd), resulting in P02 ~ 0.5. In this region the 
pulse amplitude exceeds the (0,2) singlet-triplet splitting 
energy Est, and the pulse offset is such that the (0,2) 
triplet states have lower energy than the (1,1) triplets. 
From the size of the blockaded region, and the conversion 
from detuning voltage to detuning energy e (Ae = 
AVg • 0.0676 eV/V, see supplemental material below for 
details), we find Est = 124 ± 4/ieV. 

Figure [2][a) and (b) show single-shot initialization and 
readout of (1,1) singlet and triplet states using real-time 
measurement of the charge state while pulsing across the 
(l,l)-(0,2) transition. The system is initialized by start- 
ing from the ground state S02 at < e < Est- The oc- 
cupation of S02 is verified by measuring the charge state: 
S02 is the only (0,2) state accessible at this detuning. We 
then pulse to e < to transfer the prepared S02 to the 
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FIG. 2: Single-shot initialization and readout of singlet and 
triplet states. (a),(b) Real-time measurements of Iqpc as 
the system is initialized to Sn then read out 1.7 ms later. 
We identify the final state in (a) as one of the (1,1) triplets 
(Til) because the (1,1) charge state survives for over 1ms 
during the readout. In (b) a singlet is identified because the 
system tunnels quickly back to (0,2) during the readout, (c) 
Schematic stability diagram. The points marked are the four 
detuning values used in the measurements. At _B|| > 0, Est 
is decreased by g^sB^^. The pulse is offset to keep the circle 
inside the blockaded region without changing the separation 
of the circle and triangle points. Dashed triangles bound the 
region where (l,l)-(0,2) transitions occur primarily by inter- 
dot tunneling, (d)-(g) Pulses repeatedly switch the ground 
state between (1,1) and (0,2) at 300 Hz. In (d)-(f) the system 
is often blockaded in a (1,1) triplet. With increasing mag- 
netic field from (d) to (f), the durations of blockade increase 
significantly. In (g), the pulse reaches into (0,2) far enough to 
exceed Est, and tunneling from (1,1) to (0,2) occurs freely 
for all spin states. 



(1,1) singlet Sii. To measure the (1,1) spin state at some 
later time, we pulse back to < e < Est where a sin- 
glet can tunnel quickly to (0,2) but the triplets cannot. 
The measurements are performed using detuning pulses 
with two levels that are at the positions of the filled tri- 
angle and circle in Fig.[2jc), which correspond to detun- 
ing energies of e ~ — 160 /xeV and 60 /xeV respectively at 
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We measure the lifetimes of the (1,1) singlet and triplet 
states by detecting the spin state as we repeatedly pulse 
back and forth across the (l,l)-(0,2) transition at a fre- 
quency of 300 Hz. Fig. [2](d)-(f) show real-time measure- 
ments of the charge state as the pulses are applied. In this 
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regime spin blockade is active and the system switches 
randomly between free shuttling of a singlet state and 
blockade of a (1,1) triplet state. The typical length of 
time spent in a blockaded triplet increases dramatically 
as increases. Fig[2]jg) is a control, demonstrating that 
charge shuttles freely in both directions when the pulse 
is offset to reach outside the spin-blockade regime. 

To determine the lifetimes of the states at =0 
we plot in Fig. |3ja) and (b) the number of times that 
blockaded periods of duration tb and un-blockaded peri- 
ods of duration tu are observed in 6.4 minutes of data 
(115,200 pulse periods). The histograms are very well 
fit by exponential decays, and fits to the two distribu- 
tions give lifetimes of = 9.6 ± 0.2 ms for the blockaded 
configuration and of t,j = 23 ± 3 ms for the un-blockaded 
configuration. This shows that the lifetimes of the singlet 
and triplet states are ~ 10 ms. 

The -B|| =0 lifetimes are two orders of magnitude 
longer than have been seen in comparable low-field mea- 
surements of GaAs quantum dots [H 1^ . We suggest that 
this is due to the small hyperfine coupling in natural sili- 
con, arising from the high abundance of zero-spin nuclei. 
At =0, the (1,1) triplets are degenerate and sepa- 
rated from Sii by an energy J(e) « t^/e. We expect 
singlet-triplet mixing to be driven by a small magnetic 
field difference between the two dots, resulting from the 
contact-hyperfine interaction with nuclear spins [ISHTT] . 
Predictions for the hyperfine coupling of (1,1) spin states 
are /i ~ 3 neV in silicon [T7] , compared to measured val- 
ues of /i ~ 50neV in GaAs [51 [TS]. The expected coupling 
is small enough that, in our measurements, it would be 
exceeded by the exchange splitting J. Given tc and the 
pulse amplitude, hyperfine induced singlet-triplet mixing 
should be suppressed by a factor of (1 -I- (J//i)^) ^ 500, 
compared to the maximum mixing rate when J <^ h. 

The values r„ and Tb show the time scale of singlet- 
triplet mixing, but they do not directly correspond to 
mixing times in any static configuration of the system. 
This is because the pulses continuously switch between 
two configurations, one at e < and one at e > 0. The 
singlet-triplet mixing times may be different in the two 
configurations, and at e > there are also fast, one- 
way transitions from Sn to So2- We relate the measured 
values of Tb and r„ to singlet-triplet mixing times in the 
two configurations of the system by using rate equations 
to model state occupations during a single pulse cycle. 
The inputs to the model are two times; one time r_ is 
the mixing time when the ground state is Sn during the 
e < half of the pulse, and the other time t+ is the 
mixing time when the ground state is S02 during the e > 
half of the pulse. Tunneling between Sn and S02 is 
assumed to be instantaneous. Mixing during the pulse 
transitions is ignored because the period of the pulse is 
10^ times larger than the pulse rise time. We solve for 
T+ and r_ by numerical optimization of the model to 
match the measured values of r„ and Tb (see supplemental 
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FIG. 3: Extracting lifetimes from real-time measurements, 
(a) Histogram of the number of times that the system is block- 
aded for a time tt in many measurements such as Fig. . 
The binning resolution is the pulse period. The solid line is 
an exponential fit yielding a lifetime for the blockaded config- 
uration of Tb — 9.6 ms. (b) Histogram of un-blockaded times 
(tu) for the same data as (a). An exponential fit yields a life- 
time for the un-blockaded configuration of r„ — 23 ms. (c), 
(d) Histograms of and ti, at _B|| — 250 mT. There are two 
decay rates describing blockade: at small tb there is an decay 
similar to that at zero field (r^ — 10 ms). At long tb a slower 
decay dominates (rj, = 28ms). We interpret the shorter life- 
time as that of the To, and the longer lifetime as that of the 
T_. (e) Fitted lifetimes as a function of magnetic field. The 
lifetime of blockade due to T_ states (Tb) increases with field, 
while the contribution from To and Sn states (r^ and r„) is 
field independent, (f) T_ lifetime tt, and Sh-Tq mixing rate 
at positive (negative) detuning r+ (t_). 



material below for details.) We find t_ = 24.5±3ms and 
r+ = 5.8 ± 0.3 ms. We attribute the difference between 
and T_ to a difference in tc between the two halves of 
each pulse cycle. 

As increases from T, we observe a qualitative 
change in the spin dynamics: the statistics of the block- 
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aded durations show two separate decay constants. As 
shown in Fig. |3][c) and (e), there are short blockaded pe- 
riods whose typical lifetime is field independent, and 
there are longer blockaded periods whose lifetime tj, in- 
creases with field. The two lifetimes arise because the 
system can be blockaded if it is in either a Tq or a T_ 
state, and the T_ has a field dependent energy, whereas 
the To does not. The T+ state does not play a role at 
-B|| > because its higher energy means that it is rarely 
populated. Combined with statistics of un-blockaded du- 
rations, as in Fig. |3|d), each measurement at > can 
contain simultaneously information about the lifetimes of 
three states: the Sn, T_ and the Tq. 

Fig. |3]^f) shows the T_ lifetime tt and Sh-Tq mixing 
times r+ and t_ calculated from the data in Fig.jsj^e). We 
find and r_ from t„ and using a rate equation model 
similar to the zero field case, but with no transitions to 
T+ and T_ included. This is because mixing from the 
Sii or To to the T+ and T_ will be suppressed due to 
their separation in energy. At > 0.5 T, the system 
spends so much time in the T_ state that it is impractical 
to collect enough statistics to accurately determine r^. 
Within the range of where can be measured, the 
Sii-To mixing rates are largely independent of field and 
similar to the rates seen at = 0. 

The time tt is the lifetime of the T_ during the e > 
half of the pulse and is well approximated as tt = /2 
at high magnetic field. During the e < half of the pulse, 
T_ is the ground state and it will remain populated with 
high probability when f//ii?|| > fc^T. In the e > half of 
the pulse the T_ is the first excited state and can decay 
to the So2 ground state at a of rate t^^. Such transitions 
could be induced by phonons and a spin non-conserving 
process such as hyperfine coupling O [TH [TB] or spin- 
orbit coupling [Tofe] . We find that the T_ lifetime tt 
increases strongly with field, rising to 3 seconds by = 
IT. 

In summary, we have shown that we can initialize the 
singlet-triplet qubit state into a singlet and subsequently 
measure, in single-shot mode, transitions to the (1,1) 
triplet states. Using this initialization and real-time mea- 
surement, we have measured the lifetime of singlet and 
triplet states versus magnetic field. At zero magnetic 
field, the lifetime for the singlet and all three triplets is 
~ 10 ms. At non-zero field, the To and Sn lifetimes are 
almost unchanged, whereas the T_ lifetime grows signif- 
icantly, reaching 3 seconds at IT. 
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SUPPLEMENTAL INFORMATION 
Calibration of detuning energy 

We find the conversion between detuning voltage 
and detuning energy e by measuring the width of the 
(l,l)-(0,2) transition as a function of fridge temperature 
T fridge- The relevant temperature scale for the double- 
dot system is the electron temperature Tg, which we 
model as 



fridge 



(1) 



where To is the base electron temperature of our measure- 
ment setup. In the limit where kgT^ is greater than the 
inter-dot tunnel coupling tc, the width of the (l,l)-(0,2) 
transition will be determined by Te- When fcsTe <^ tc, 
the width is independent of temperature (TUl . 

The line-shape of the inter-dot transition was measured 
by detecting the charge-sensor current Igpc with a lock- 
in amplifier while gates 2 and 4 were swept across the 
transition in the direction (as shown in Fig. 1(b) of the 
main text): AV2 = — 1.375AV4. The lock-in excitation 
was also applied in the Ve direction by modulating gates 
2 and 4 with a relative phase of 180 degrees and relative 
amplitudes of V^p — 0.35V2p. The excitation frequency 
was 500 Hz and the amplitude of the excitation at the 
gates was ^ 20 /iV on gate 4. The results were fit with 
the expected line-shape in the thermally broadened limit, 
which is the derivative of a Fermi function: 



= Acosh-2 



ajVe - Veo) 

2kBTc 



(2) 



where the voltage changes along the detuning axis are 
given by = V'e and V2 = — V'e/1.375, and a is the 
conversion factor between and detuning energy. 



Figure SI shows the measured values of the transi- 
tion width kBTc/a as a function of T fridge- We find that 
the width depends on temperature down 50 mK, confirm- 
ing our assumption that we are in the limit fc^Te > t^, 
and also indicating that the base electron temperature 
To is < 50 mK. By fitting Equation [l] to the data, 
we find the calibration for the detuning energy to be 
a = 67.6meV/V. 



Zero magnetic field rate equation model 

A set of coupled rate equations are used to extract Sn- 
Tq mixing rates at zero magnetic field from the measured 



5 




0.2 

Tfridge (K) 



FIG. SI: Measured widths of the (l,l)-(0,2) transition as a 
function of the fridge temperature Tfridge- The soHd line is 
a fit of the form given by Equation [l] with To ~ 10.3 mK. 
The slope of the dependence at large temperatures yields the 
conversion between gate voltage along the detuning axis and 
energy: a = 67.6meV/V. The dotted line shows the con- 
version if Tq is neglected, which agrees well with the data at 
large temperatures. 



values r„ and Th, which are the hfetimes of blockaded and 
un-blockaded configuration of the system. We relate 
and Th to the results of the rate equation model via two 
values: Pr = e~^l'^^ is the probability of remaining un- 
blockaded in a single pulse cycle (the probability of re- 
turning to (0,2)), and P/v = er'^l'^^ is the probability of 
remaining blockaded after a single pulse cycle (the prob- 
ability of not returning to (0,2)). T is the period of the 
pulse cycle, which in all the experiments consists of two 
halves of equal duration. For the first half of the cycle 
(t < r/2), the system is at negative detuning (e < 0), 
which favors the (1,1) state, while during the second half 
of the cycle (T/2 < t < T), the system is at positive 
detuning, where the (0,2) state is favored. 

The physical mechanisms we consider are inter-dot 
tunneling between Sn and S02, and singlet-triplet mixing 
due to a small field difference between the dots arising 
from the hyperfine interaction with background nuclear 
spins. At zero external field, the time-averaged effect of 
the small field difference between the dots is to cause all 
the (1,1) triplet states to mix with the (1,1) singlet at the 
same rate [T21 [H] . This mixing rate will depend on the 
singlet-triplet splitting J « i^/e, where tc is the inter-dot 
tunnel coupling and e is the detuning energy. The two 
halves of each pulse cycle will have different values of |e| 
and, possibly, slightly different values of tc- We there- 
fore allow for two singlet-triplet mixing rates for the two 
halves of the cycle. 

The double-dot states included in the model are Sn, 
T_, To, T+ and S02, with occupation probabilities pi, 
P2, P3, P4 and P5 respectively. During the first half of 
the pulse cycle, while the system is at negative detuning 
(e < 0), we allow all three triplet states to mix with Sn at 
a rate p_ . The S02 state is assumed to be un-populated 



because it will quickly tunnel to Sn. The evolution of 
the occupation probabilities is given by: 



(3) 
(4) 
(5) 
(6) 
(7) 



During the second half of the pulse, from t = r/2 to 
t = T, the system is at positive detuning (e > 0) and 
the S02 is now the ground state. The rate of singlet- 
triplet mixing is and we assume that all occupation 
of Sii is instantly transferred to S02. The evolution of 
the occupation probabilities is given by: 
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We note that the model yields the same results if tran- 
sitions between triplets are also allowed. 

The two sets of coupled rate equations, Eqn. [3][7] and 
Eqn. [8p2l are solved numerically to predict Pn and Pu 
for various p- and p+. The values of p_ and p+ are 
then optimized to best match P^ and Pn given by the 
measurements. Errors on p_ and p+ are found by re- 
calculating the model multiple times while varying Pj^ 
and Pr within a range given by the errors on and t,i. 

To calculate Pr, the initial conditions are set to unit 
probability of being in Sn: Pi(0) = 1 and pjv^i(O) = 0. 
The occupation probabilities are evolved numerically for 
a time T/2 according to Eqn. 3][7 The resulting values 
of pn are then used as initial conditions for a further 
evolution for time T/2 according to Eqn. 8p2 however, 
all occupation of Sn is first transferred to S02 to simulate 
the fast inter-dot tunnel coupling of these states. The 
probability of remaining un-blockaded for a single pulse 
period is then given by Pr — (T) . 

Pn is calculated in a similar way to Pr, but with the 
initial conditions being a unit probability of starting in a 
triplet state. The probability of remaining in a blockaded 
state for a single period is given by P/v = X]n=2P"(-^)- 
Because of the symmetry of the rate equations, the an- 
swer is independent of which triplet is initially occupied. 

The lifetimes found by experiment at zero magnetic 
field are r„ = 23±3ms (Pr = 0.87) and n = 9.6± 0.2 ms 
(Ptv — 0.71). We find that these values are consistent 
with Sii-To mixing times of 1//?+ = 5.8 ± 0.3 ms and 
l/p_ = 24.5 ± 3 ms. 
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